Generalized Master Function Approach to 
Quasi-Exactly Solvable Models 

M. A. Jafarizadeh a ' 6 ' c *, S. J. Akhtarshenas a ' b ' c f 

Department of Theoretical Physics and Astrophysics, Tabriz University, Tabriz 51664, Iran. 
^Institute for Studies in Theoretical Physics and Mathematics, Tehran 19395-1795, Iran. 
c Pure and Applied Science Research Center, Tabriz 51664, Iran. 

February 7, 2008 



*E-mail:jafarzadch(Q!ark.tabrizu.ac.ir 
tE-mail:akhtarash@ark.tabrizu.ac.ir 



1 



Master Function Approach to Quasi-Exact Models 



2 



Abstract 

By introducing the generalized master function of order up to four together with corresponding 
weight function, we have obtained all quasi-exactly solvable second order differential equations. It is 
shown that these differntial equations have solutions of polynomial type with factorziation properties, 
that is polynomial solutions P m (E) can be factorized in terms of polynomial P n +i{E) for m > n+ 1. All 
known quasi-exactly quantum solvable models can be obtained from these differential equations, where 
roots of polynomial P n+ i(E) are corresponding eigen- values. 

Keywords: Quasi-Exact, Differential Equations, Orthogonal Polynomials, Factoriza- 
tion. 

PACs Index: 03.65.Ge 
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1 INTRODUCTION 

During the last decade a remarkable new class of quasi-exactly solvable spectral problems was introduced 
0, ||, ||, [|, II . These occupy an intermediate position between exactly solvable and unsolvable models in the 
sense that exact solution in an algebraized form exists only for a part of the spectrum. 

The usual approach to the analysis of quasi-exactly solvable systems is an algebraic one in which the 
operator is expressed as a non-linear combination of generators of a Lie algebra. Another recent developemcnt 
is the work of Bender-Dunne where they have shown that the eigen- functions of a quasi-exactly solvable 
schrodinger equation is the generating function for a set of orthogonal polynomials P m (E) in energy variable. 
It was further shown that, these polynomials satisfy the three-term recursion relation. Also, all polynomials 
beyond a critical polynomial P m {E) factorize into the product of polynomial P n+ \{E) and another arbitrary 
polynomial. 

In this paper we suggest a generalization of Bender-Dunne approach to all possible one-dimentional 
quasi-exactly second order differential equations. 

For this purpose, the succesful master function approach of references (?], || to exactly solvable models, 
is generalized to a master function of order up to four which gives all possible one-dimensional quasi-exactly 
solvable models, where Bender-Dunne model || and Heun differential equation H are among them. 

The paper is organized as follows: In section II we show that we can generalize the usual quadratic 
master function to a master function of at most four order polynomials, then the most general quasi-exactly 
solvable differential operators related to generalized master function of degree k = 3 and k — 4 are given, 
respectively. 

In section III, expanding their solutions in powers of x, we get 3-term and 4-term recursion relations among 
their coefficeints, where Bender-Dunne factorization follows through imposing the quasi-exactly solvability 
conditions of section II. At the end of this section we list all possible related quasi-exactly solvable differential 
equations for k — 3 and k — 4 in Tables I and II, respectively. 

Finally at section IV, we derive all possible one-dimensional quasi-exactly solvable quantum Hamiltonian 
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from the differential operators of section III, via prescription of references 0, where we have listed them 
at the end of section III, except for those which can given in terms of elliptic functions. Paper ends with a 
brief conclusion. 



2 QUASI-EXACTLY SOLVABLE DIFFERENTIAL EQUATIONS 
ASSOCIATED WITH GENERALIZED MASTER FUNCTIONS 

By generalizing master function of order up to two J?], |8| to polynomial of order up to k, together with 
the non-negative weight function W(x), defined at interval (a,b) such that w ^ x ^ (A(x)W(x)) to be a 
polynomial of degree at most (k — 1), we can define the operator 

where B(x) is a polynomial of order up to (k — 2). The interval (a, b) is chosen so that, we have A(a)W(a) = 
A(b)W(b) = 0. 

It is straightforward to show that the above defined operator L is a self adjoint linear operator which at 
most maps a given polynomial of order m to another polynomial of order (m+k — 2). Now, by an appropriate 
choice of B(x) and weight function W{x), the operator L can have an invariant subspace of polynomials of 
order up to n. Then by choosing the set of orthogonal polynomials {4>o(x), 4>\{x), • • • , (f> n (x)} defined in the 
interval (a, 6) with respect to the weight function W(x): 

b 

(f> m {x)(f>n{x)W(x)dx — 0, for m — n (2-2) 
as the base, the matrix elements of the operator L on this base will have the following block diagonal form: 

Lij =0, if {i < n and j ' > n + 1} or {i > n + 1 and j < n}. (2-3) 
Since, according to the well known theorem of orthogonal polynomials, 4> n {x) is orthogonal to any poly- 
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L = 



M 
N 

where M is an (n + 1) x (n + 1) matrix with matrix elements 



(2-4) 



Mij= dxW(x)(j) l (x)L(x)(j) j (x), i,j = 0,1,2, ...,n, (2-5) 

</ a 

and JV is an infinite matrix element defined as above with i, j > n+ 1. 

The block diagonal form of the operator L indicates that by diagonalizing the (n + 1) x (n + 1) matrix 
M, we can find (n + 1) eigen-values of the operator L together with the related eigen-functions as linear 
functions of orthogonal polynomials {(f>o(x), (f>i(x), ■ ■ ■ , 4> n {x)} ■ 

In order to determine the appropriate B(x) and W(x) for a given generalized master function A(x), we 
Taylor expand those functions: 



A«(o) 



i=0 



where A® (0) = \ x=0 



(A(x)VK(x)) _ g ( — 



(0) 



W(x) 



where 



i=0 



{AWy_\ W = d { W( X ) ) 
V W ) 1 ' dx l 



\x=0 



k 2 B^(0) 
B(x) = Y^ r^x\ where 5 W (0) 



i=0 



d l B(x) 
dx i 



\x=0 ■ 



(2-6) 



(2-7) 



(2-8) 



Then, the existence of invariant subspace of the polynomials of order n of the operator L leads to the 
following linear equation between the coefficients of above Taylor expansions: 



^(»+2) 



1(1-1) 



(AW)' 
W 



(i+1) 



(i + iy. 



B(0 
^ + — r- = 0, 



(2-9) 
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where 

I = n, and i = 1, 2, fc — 2 

Z = n — 1, and i = 2, 3, fc — 2 

Z = n — k + 4, and i = k — 3, fc — 2 
Z = n — fc + 3, and i = k — 2 

The number of above equations, for a given value of fc, is ( fc-1 )( fc ~ 2 ) _ jf we arc to determine only the 
unknown function B(x) without having any further constraint on the weight function W(x), then the above 
(fc-i)(fc-2) e q ua tj ons should be satisfied with (fc — 2) coefficients of Taylor expansion of B(x) as the only 
unknowns, since can be absorbed in the eigen-spectrum operator L. Therefore, we are left with (fc — 2) 
unkowns to be determined, where the compatibility of equations (2-9) require fc — 3 at most. On the other 
hand, if we add the coefficients of Taylor expansions of A(x) and ^^^^^ to our list of unknowns, ( to 
be determined by solving equations (2-9) ), then their compatibility conditions require that: 

(fc-l)(fc-2) 



3(fc - 1) > 



(2-10) 



or fc < 8, where further investigations show that we can have at most fc = 4, since for fc > 5 the coefficients 
A^(0) and ( ^^p ^ (0) will vanish. Below we summarize the above-mentioned discussion for fc = 3 
and fc = 4, separately. 



2.1 A; = 3 

In this case, B{x) is a second order polynomial where B^ can be determined by solving equation (2-9): 



(2) 



(2-11) 



which is the only unknown in this case. 



Master Function Approach to Quasi-Exact Models 7 
2.2 k = 4 

Again, the solving of equations (2-9) leads to: 

4(4) 

B^=-—n{n-l), (2-13) 

and 

Here, besides having constraint over second order polynomial B(x), we have to put further constraints on 
the weight function W(x) given in (2-14). 

Dcfinctly we can determine n + 1 eigen-spectrum of the operator L, simply by diagonalizing the [n + 
1) x (n + 1) matrix M, since it is a self-adjoint operator in Hilbert space of polynomials and it has a block 
diagonal form given in (2-4). 

As we are going to see in the next section, we can determine its eigen-spectrum analytically using some 
recursion relations. 



3 RECURSION RELATIONS 

In this section we show that the eigen-functions of the operator L are a generating function for a new set 
of polynomials P m (E) where the eigen- function equation of the operator L leads to the recursion relation 
between these polynomials. Quasi-exact solvable constraints (2-9) will lead to their factorization, that is, 
P n +N+i{E) = P n +i(E)QN{E) for N > 0, where roots of polynomials P n+ \{E) turn out to be the eigen- values 
of the operator L. 

To achieve these results, first we expand 4>(x), the eigen-function of L, as: 

oc 

1>(x) - ]T P m {E)x m , (3-15) 

m=0 
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where eigen-function equation: 



Lip(x) = Exp{x) 



can be expressed as: 



-A(x)J2m{m-l)P m (E) 



„ 2 _ (A {x )w {x) y g mPm[E)xm -, 



W{x) 



m—1 



+ B( X ) J2 Pm(E)x m = EJ2 P™{E)l 



m— m— 

which leads to the following recursion relations for the coefficients P m (E): 

x'\ (°) 



A™ (m + l)(m + 2) + (^^-) (m + 2)^ P ro+2 (P) 

i)+(^^) () (m + i) + ^p m+1 (i;) 



MP) 

+ — m(m+l) 



P m (E) 



) 



+ 



/\(3) 



V 



— (m - l)(m - 2) + A ^—m - — | P ni . ,.(£)=(). 



Below we investigate recursion relations thus obtained for k = 3 and k — 4, separately. 



3.1 A; = 3 



In this case the 4-term general recursion relation reduse to the following 3-term recursion relation: 

I \ (o) 



A« ( m + l)(m + 2) + (^p-^j (m + 2)^ P ro+2 (P) 



A(2) 



+ [ ^ r m(m+l)+ (^^) ( \m + l)+p)p m+1 (P) 



——m(m - 1) + ± — J- 



+ | ^mlm - 1) + v " 2! / m - B« | P m (E) = 0. 
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In order to have finite eigen-spectrum, that is, quasi-integrable differential equation, the above recursion 
relation should be truncated for some value of m = n, which is obviously possible by an appropriate choice 
of: 

n ( 4(3) f(AW\'\ (2) \ 

(3-20) 



which is in agreement with the result of previous section given in (2-11). 

Using the recursion relation (3-19), with given in (3-20), we get a factorization of polynomial 

Pn+N+i{E) for N > in terms of P n +i(E) as follows: 

Pn + N+i(E) = P n+1 (E)Q N (E) N>0 (3-21) 

where, by choosing the eigen- values E as roots of polynomial P n +i(E), all polynomials of order higher than 
n will vanish. 

In order to determine corresponding eigen- functions, it is sufficient to evaluate P m (Ei) for m = 0, 1, 2 ... n 
with Ei as roots of P n+ i(E), then eigen- function ipi(x) corresponding to eigen- value E t can be written as: 



ij>i{x) = J2 P m {E l )x m , t = , 1 ... ,n. (3-22) 



The above eigen- functions are polynomials of order n, hence they can have at most n roots in the interval 
(a, 6), where, according to the well known oscillation and comparison theorem of second-order linear differ- 
ential equation Jl0| , these numbers order the eigen- values according to the number of roots of corresponding 
eigen-functions. Therefore, we can say that the eigen-values thus obtained are the first n + 1 eigen- values of 
the operator L. 

Using the recursion relations (3-19), we can evaluate the polynomials P m (E) in terms of Pq(E), where 
we have chosen Pq(E) = 1. We have evaluated the first five polynomials which appear in the Appendix I. 

As an illustration we give the results for A(x) = x and n = 3 with a = 1, (3 = 0, 7 = —1 which is 
equivalent to the Bender-Dunne model: 

Px{E) = -\e, 
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P 2 {E) = -l + ^E\ 



p ( E ) = -E-—E 3 , 
iy '4 144 ' 



PJE) = — - —E 2 + —E 4 . 
4V ' 10 48 2880 

Obviously P m (E) have the parity of m. 

By finding the 4-roots of Pa(E) we determine the corresponding four eigen- values: 

E = -7.398556194, 

Pi = -2.293766823, 
E 2 = 2.293766823, 
E 3 = 7.398556194. 

Finally for the coefficient P m (Ei) we get: 



P (E ) = 1, Pi (E ) = 3.699278097, P 2 {E ) = 3.561552813, P 3 (E ) = -.962769686, 
P {) {Ei) = 1, Py(Ex) = 1.146883412, P 2 (E 1 ) = -.5615528135, P 3 (^i) = -.4896337383, 
P (-E , 2 ) = l, Pi(J5 2 ) = -1.146883412, P 2 (S 2 ) = -.5615528135, P 3 (£ 2 ) = .4896337383, 
P (J5 3 ) = 1, Pi (£? 3 ) = -3.699278097, P 2 (P 3 ) = 3.561552813, P 3 (P 3 ) = -.962769686. 

Using the above coefficients we can determine the corresponding eigen-functions through formula (3-22) . 

In Table I we give all quasi-exactly solvable operators which can be obtained by choosing different 
generalized master function of order 3. This Table contains all possible models corresponding to different 
choice of A(x) up to translation and rescaling of variable x. Also by choosing A[x) as a polynomial of up to 
second order with 7 = we lead to the exactly solvable models of references R pf . 
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A(x) 


W{x) 


L(x) 


intervals 


X 


< x < +oo 
a > -1, -oo < (3 < +oo, 7 < 


-x£? -(a + l + (3x 
+2-/x 2 )£ + 2n-fx 


x 2 


x a> e 0/x+jx 

< X < +oo 
— oo < a < +oo, (3 < 0, 7 < 


-x 2 £, + ((3-(a + 2)x 
-7^ 2 )s + n ^ x 


x(l — x) 


x a (l -xfe-^ 
< x < +1 
a > — 1, /3 > -1, -oo < 7 < +oo 


x{x + {-a -\ 
+ {a + [3 + j + 2)x 
-lx 2 )-^ + n^x 


X 3 


x a e -l3/x 2 -~f/x 

< x < +oo 
a < -3, /3 > 0, -oo < 7 < +oo 


- x z£-(2f3 + 1 x 
+(a + 3)a; 2 )£ + n(n + a + 2)x 


x 2 {l-x) 


x a {l -xfe-~<l x 
< x < 1 
-oo < a < +oo, f3 > -1, 7 > 


x 2 (x-l)^ + (-7+( 7 - a -2)x 
+ {a + (3 + i)x 2 )± 
-n(n + a + [3 + 2)x 


x(l + x 2 ) 


x a (l + x 2 ) /3 e 7 * a "' la; 
< a; < +oo 
a > -1, f3 < -(a + 3)/2, -oo < 7 < +00 


-x(l + x 2 )£s - (a + l+jx 
+ (a + 2f3 + 3)x 2 )£ 
+n(n + a + 2(3 + 2)x 


x(l — x)(a — x) 
a > 1 


^(l-z^a-a;) 7 
< x < +1 
a > -1, (3 > -1, -00 < 7 < +00 


x(x - l)(a - x)£, + {-a(a + 1) 
+((a + l)(a + 2) + a/3 + j)x - (a + (3 
+7 + 3)x 2 )£ + n(n + a + (3 + 7 + 2)x 



Tabic 1: Quasi-exactly differential operators obtained from generalized master function of order up to 3 
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3.2 k = A 

Again in order to truncate the recursion relation (3-18) and to factorize polynomials P n+ N + \(E) in terms 
of P n+ i(E), we should have: 



n 



2 \ 3 



AGO _ AW, u , (^) 



(3) 



2! 4! v M ' 3! 

and 



(n-l)(n-2) + -* ^^n, (3-24) 



(3) 

B (2) A< 4 ) 



(AW)' 
W 



2! 4! ' " 3! 
Solving the above equations we get: 



n(n-l) + - ^ — (n + 1). (3-25) 



4( 4 ) 

S (2) =-— n(n-l), (3-26) 



and 



/ (Awy \ (3) .i Mi 



(n-l). (3-27) 



V W y 2 

The equations (3-23) , (3-26) and (3-27) are the same equations which are required in the reduction of 
the operator L to its block diagonal form. 

Again roots of polynomials P n +i(E) will correspond to n+ 1 eigen- values of the differential operator L 
with eigen-functions which can be expressed in terms of P m (Ei) for m < n, where polynomials P m {E) can 
be obtained from recursion relation by choosing Pq(E) = 1 and P-i(E) = 0, where we have given the first 
four polynomials in Appendix II. 

In Table II we list all quasi-exactly differential operators which can be obtained from the generalized 
master function of order up to four. 

Tables I and II contain all quasi-exactly second order differential equations which can be obtained from 
Lie algebraic methods Q. For example we get the Bender-Dunne model || for the choice of A(x) = x and 
(3 = 0, and similarly, we get the Hcun differential operator ( Fuxian equation with four regular singular point 
) ||, for the choice of A(x) = x(x — l)(x — a). 
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A(x) 


W(x) 


L(x) 


intervals 


x 4 


x a e l3/x 3 +j/x 2 +S/x 

< x < +oo 
(3 < 0, — oo < 7 < +oo, — oo < 5 < +oo 
a = -2(n+ 1) 


+<5x 2 - (a + 4)x 3 )^ - nSx - n(n - l)x 2 


x 3 (l - x) 


x a (l — x)@e~ J / x ~ s / x2 
< x < +1 
(3 > -1, -oo < 7 < +oo, (5 > 
a = -2(n + l) - /3 


x3( x -1)4 + (-25 + (25 - y)x 
+ (7 - a - 3)x 2 + (a + 13 + 4)x 3 )^| 
+n(n + a — 7 + 2)x + n(n — l)x 2 


x 2 (l + x 2 ) 


_i_ t,2\/3 7/x+<5tan _1 (a:) 

< x < +oo 
— oo < (3 < +oo, 7 < 0, — oo < 5 < +oo 
a = -2(n + (3+l) 


-x 2 Cl + x 2 )-£y + h - (a + 2)x 
+ ( 7 -<5)x 2 -(a + 2/? + 4)x 3 )£ 
+n(5 — 7)x — n(n — l)x 2 


T f 1 — T\(o — t) 

a > 1 


x a (l - x^fa - xVe 6 ^ 
< x < +1 
/3 > -1, -oo < 7 < +oo, 5 < 
a= -2(n + 1) - (3 - 7 


r 2 (r — l)(a — T~)-£r 4- (a,5 — (an> 4- (a 4- 1)5 

+2a)x - ((a + l)a + af3 + 7 + (5 
+3(a + l))x 2 - (a + 13 + 7 + 4)x 3 )^ + n((a 
-n + 4) (a + 1) + a/3 + 7 + (5)x - n(n - l)x 2 


x 2 (l - x) 2 


x a (l - x) p e l/x+s ^ 1 ^ x) 
< x < +1 
-00 < 13 < +00, 7 < 0, (5 < 
a = -2(n + l) - (3 


-x 2 fl - x) 2 -£? + h- - (a + 2i + 2)x 
+ (2a + (3 + 7 - 5 + 6)x 2 - (a + (3 + 4)x 3 )£ 
-n(+2n + 2a + (3 + 7 - 5 + 4)x - n(n - l)x 2 


x(a - x)(l + x 2 ) 
a > 


x a {a - x)? (1 + x 2 )"< e stan ~ lx 
< x < a 
—00 < 7 < +00, —00 < 5 < +00 
-l</3< -In- 2 7 - 1 
a = -2(ra + 1) - /? - 27 


x(x - a)(l + x 2 )^ + (-a(a + 1) 
+ fa + /? - aC5 - 2))x - (a(a + 2-7 + 3) 
-5)x 2 + (a + /3 + 2 7 + 4)x 3 )£ 
+n(a(n + a + 27 + 2) - <5)x + n(n - l)x 2 


x(l — x)(a — x)(b — x) 
1 <a<b 


x a (l-x)^ f (a-x)T(6-x)' 5 
< x < +1 
— 1 ^ R <r On — -v — 1 


x(x- \){a - x)(b - x)£t 
+ (-ab(a + 1) + (2(o + b + ab) + a6(a + /3) 

4- -v"\ -I- /if'rv 4- AWr ((n 4- ?i Vrv 4- R 4- 



Master Function Approach to Quasi-Exact Models 14 

4 QUASI-EXACTLY POTENTIAL ASSOCIATED WITH GEN- 
ERALIZED MASTER FUNCTION 

As in references ]7[ ||, writting: 

iP(t)=A 1 / 4 {x)W 1/2 (x)(f>(x), (4-28) 

with a change of variable jj? = y/A(x), the eigen- value equation of the operator L reduces to the schrodinger 
equation: 

H(t)1>{t) = Ej>(t), (4-29) 

with the same eigen- value E and i/j(t) given in (4-29), in terms of the eigen function of L, where H(t) = 
' V(t) is the similarity transformation of L(x) defined as: 



dl 



with: 



H(t) = A 1 ' 4 {x)W 1/2 {x)L{x)A- 1 ' 4 {x)W- 1/2 {x) (4-30) 



vm 3 A3 ® _ + i Mt)w(t) , \Mt) i w(t) 

[> ~ lGA^t) 4W 2 (t) + 4A(t)W(t) + 4A(t) + 2W(t) + [ >- [ ' 



It is also straightforward to show that: 

dt<t>(t)H(t)ip(t) = ( dxW(x)ip(x)L(x)ip(x), (4-32) 

J a 

Hence block diagonalization of L leads to block-diagonalization of H . 

As an illustration we give below an example with A(x) = x 3 , weight function W(x) = x a e 13 ] > lx ' 1—1 '' >x and 
interval [0, oo), where a < —3 and (3, 7 > 0. 



From a change of variable dx/dt — -JA(x), we get x(t) = A/t 2 , hence using equation (4-31) we have for 
potential V(t): 

1 , ^ f 15 9 , 9 „ \ 1 1 / „ 1 ,\ , 



= 7; {a + I) + {— + a 2 + Ana + Aa + An 2 + 8n j — + - la f3 + ~Y ) t 



Master Function Approach to Quasi-Exact Models 



15 



Below we give a list of quasi-exact solvable potentials, except for those potentials which can be expressed 
in terms of elliptic functions, since in this case we get rather long expressions for them: 



A(x) = x, x(t) = — 



V n (t)= (a + l)+ [a 2 - 1 



2 v ' \ A t 2 



A(x)=x 2 1 x(t)=e t 



V n (t) = i (1 + a 2 - 2/3 7 + 2a - 2ape~ t + f3 2 e- 2t + 7 2 e 2t + 2 (2 7 + 2n 7 + a 7) e 4 ) , 



4 / x /, \ / \ 1 + sin(t) 
A(x)=x(l-x), x(t) = — ^ 



K(*) = ^+n 7 -a/3- / 8-a+i (/? 7 -a 2 - /3 2 -a 7 -l) + ^ + 7 + ^ + sin(t) 



A{x)=x 3 , x(t) = ^ 
V n (t) = l(a + l)+ ^+a 2 +4na + 4a + 4n 2 + 8n^ + i /? + 2_ ) / 

+ 16 + 64* ' 



A(x) = a; 2 (l - x), x{t) = I - tank 2 {-) 

VJI) = ( ^/^ft) _ 1 ) (^-(^2n 2 + 2 + 2na+y+4n + a/3 + 2a+^ + 2n/3 + 2/3S 

1 /-7 2 «7 1 . a 2 \ , 9 a 7 , ^ 7 2 ,4 

+ - — ■! ^ + 7:-7+ a -/37+Tr cos/l + —rcosh 6 + j-cosh 4 



2 4 22' ''27 4 16 
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A(x) = x 4 , x(t) = j 
V n (t)=(^+ 1 + 2n 1 )+( 1 5 + 3nP + 3f3)t+(%^+ 1 2 ) t 2 +3(3 7 t 3 + ^t\ 



4 



A{x) = x 3 {\ - x), x(t)= 1 



4 + t 2 



V n {t) = -(-L+5 + 05+ Y+n 7 ) + (/3 2 4 , /2 



= x 2 (l + x 2 ), a;(i) = — ' 



sinh{t) 



r„(/) = |n + n 2 -^ + i + 2n/3 + /3 + /3 2 -(n7 + 7 + /3 7 )sm/i(t) 



A(ar) = :e 2 (1-:e 2 ), a;(t) ' 



cosh{t) 



..... , 1 7 2 * 2 /37 7* 1 ,^/5 2 

-(_+n/3- T + n 7 +- + - + - + — - y +n + T + n + y j «*fc(t) 
+ ( ^ + coafc(t) +(s + nS+^ + ^f] cosh 3 (t) + ^cosh 4 {t) 



2 2 2 2 / w V 2 2/ w 4 



^(x) =x 2 (l-x) 2 , a;(i) = 



1 + e* 
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-3/3<5+^-3n7 + 7 + 4n + (5 2 - 2 7 (5 + /3 2 + l + 4n 2 + 2 > 4//J | S'- 
+ (3/34-6n+^-37 ( 54-| + ^4-^4-6n 2 -2/35 + 47 4-6n/3 + 2/37-2n 7 )5' 8 



2 '244 
+ (4n 2 - 4^-27(5 + 3/37 + 1 + 2 n 7 + 4n + /3 2 + 2/3 + 7 2 + 6 7 + 4n/3 I S* 9 



'—7(5 3t 2 1 /3 2 /3 

+ ( ~y~ + ^-+47 + n + 2/3 7 + 3n7+- + ^- + ^+n/3 + n 2 )S'io 

+ f n 7 + 7 2 + 4^ + j) S u +^rSi2 



where Sk is defined as: 

S k = exp(-2 ^ + ( 7 e-'+7 + «S(e t + l))^ , fc = 4,5, • ■ ■ , 12. 

5 CONCLUSION 

As we saw by introducing of master function A{x) as a polynomial of order at most four, we could obtain all 
quasi-exactly second order differential equations. It is shown that the eigen-equation relation L^l(x) — E^>(x) 
generates a set of polynomials P m (E) , where these polynomials satisfy 3-term and 4-term recursion relations 
for master function of at most three and four, respectively. Finally the quasi-exactly solvablity leads to 
factorization of polynomials P n+ N+i{E) for N > in terms of P n+ \(E), where by determining the roots of 
P n+ i{E) we can determine first n + 1 eigen- values of these quasi-exactly solvable differential equations. 



APPENDIX I 

The First Five Polynomials P n (E), For k = 3 
To abbreviate, we set F« = (^) W , 
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E 

Pi(E) = - " 



_ 1 EF^I + E 2 + B^ 



P 3 (.E) = (-A^ EF^ - A™ e 2 - A™ F W _2EFW 2 -3F^E 2 -2 F^ B^ F^ 
-3 EB™ FW + EFW A^ + EFW -2EB^ A*')) 
/(6F(°) (aW+F<V) (2AW+FM)) , 



E A 



P 4 (E) = - (-3 A^ 2 E 2 - 3 B« 2 F<°> 2 - 6 B^ 2 F<°) A™ + 3 F™ F<°) 2 + A^ F<°) 2 

+j4 (S) S 2 F (0) + 2 A (S) ^(i) _ g E 2 B (1) A (l) +AE 2p{2) p{0) + 7 S 2 F (2) _ g £ 2 B (l) p(0) 

_ 4 S 3 _ g £ F (i) 3 _ n V _ g p(l) E 3 _ QA (2) p(l) B {1) p{0) _Q A (2) EB (1) A (l) 

+3A (2) E p(2) p(0) +iA (2) E p{2) A (l) _ 1QA (2) EB (1) p(0) _ 12 p(l) EB (1) A (l) + Qp(l) E p(2) p(0) 

+9F^ EFW A^ -UF^ EBW F^+2A^ B ^ F^ A^ + A^ EF^ F^ 
+2 A^ EFW A^+6 F^ B^ A^ - £ 4 ) 
/(mfW (aW+fW) (2AW + FW) (3^ J ) + F(°))) , 
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P 5 (E) = - (-46EAW FM A^ - 88EF^ A^ + 16 A^ EF™ A^ F^ 

+2AEF^ 2 A<-^ - 48 F^ EB^ A^ 2 + 18 EF^A^ 2 + 6 V°> 2 + 24EB^ 2 A^ 2 
+ 12 A^ EFW A^ 2 +4 AW EF^ F^ 2 - 24 A^ EB^ A^ 2 - 32 F^ A^ B^ F<°> A^ 
-60 F^ F< 2 > B< J > F<°> A^ + 50 EB^F^ A< J > - 25 EF^ B^ F^ 2 - 13 EA^ B^ F<- ^ 
-46 F^ A^ E 2 A^ + 80 F< J > E 2 B^ A^ - 40 E 2 F^ F^ - 91 E 2 F^ A^ 
+50 F^ E 2 B^ F^ _ 54,4(2) EF™ F (0) - 84 A™ F^ EF™ A^ + 137 A™ F< J > EB^ F^ 

-24 A^ A<*> B^ F<°> A^ - 10 A^ A™ EF™ F^ - 24 A^ A^ EF^ A^ 
-54 A^ F^ B^ AW + 20 F^ flW V°> 2 - 5 A<*> E 3 F^ - 14 A<*> E 3 A^ + 20 E 3 B^ A^ 
-10E 3 F^ FW - 25F 3 F^ A^ + 10F 3 F^ F^ + 15 EB^F^ 2 + 66 A<*> E 2 B^ A^ 
-33 A^ E 2 F^ F^ - 63 A<*> E 2 F^ A^ + 50 A^ E 2 B^ + 72 A<*> F^B^ F^ 
+72 FW 2 EBW A^ - 36 FW 2 EFW F (0) - 72 F^EF^ A<*> + 70 F^EB^ F^ 
-12 A^ EFW 2 F<°) -32A< s > EFW 2 AW + 48 F^ 1 ^ F^V ' A< j > - 24 F^ F^ flW F^ 2 
-12F (J) A (3) F (i) F (0)2 - 17 F (J) ^ 3 )f 2 F (0) - 24A (2) F {2) B {1) F (0)2 - 10 A^ 2) A {3) F 2 F (0) 
-24 A^ A^ E 2 A^ - 10 A^ A^ B^ F^ 2 + 108 A^ F^ EB^ A^ 
+ 10 F^E 4 + 18 A^E 2 + 27 A^^V + 10 A^ E* + 24EF ( - 1 ^ + 50 V 
+35 F( J ) 2 F 3 + E 5 + 93 A( 2 ) 2 F( J ) E 2 + 66 A( 2 ) 2 FF^) 2 + 18 A^B^ F (0 > + 18 A^EF^ 
+22 A( 2 ) S( J ) 2 F( ) 2 + 72 A^ EF^ 3 + 127 A^ F^ 2 E 2 + 65 F^ F 3 
+24F( j ) 3 F.( j )f(°)+66 j 4( 2 ) 2 f( j )f( j )f(°)+36A( 2 ) 2 FF( j )A( j ) 
-18 A( 2 ) 2 FF( 2 ) F(°) - 18 A( 2 ) 2 FF( 2 ) A( j ) + 63 A(*) 2 f;B( J ) f( ) + 48 A^ B^F^ A^ 
/(l20F(°) (a^ + fW) (2^) + F(°)) (3A( j ) + F( )) (4^>+f<»')) . 

APPENDIX II 
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The First Four Polynomials P n (E), For k = 4 



Pi(E)= 1 



l EF^ +E 2 - 
2i - 2 F(°) (At 1 ) + F(°)) ' 



P 3 (E) = (A^ EFW + A^ E 2 - A<*> B^ F^ + 2 EF^ 2 + 3F^E 2 -2 F^ B^ 
^E 3 + EB^ FW - EFW A^ - EF& F (0) + 2EBW A^ + B<*> F^ A^ + B<*> F^ 2 ) 



P 4 (E) = - (3 A^ B<«) F^ 2 + A^ B W F^ 2 + 9 A^ EF^ 2 - 3 A^ EF™ A^ 
+2A^ EBW F<°> -9A^ F^ B^ - 2EB™ F^ 2 + 3A^ B<*> F^ A< J > + 6 A<*> EB^ A™ 

_ 3 EF (2) p (0) _ 9F (D EF (2) A (D +AF (1) EB {1) p(0) _ SB (D 2 F (0) 2 _g EB (2) p{0) A (l) 

-2.4 (3) EF^ A {1) - A {3) FF (J) F (0) +2.4 (3) B {1) F {0) A {1) + 6 F (2) F {0) A (i) + 3 EF {3) A {1) F (0) 
-6 F^ EFW F (0) + 3 A<*> V + 4 A^ E 3 + 6 EF^ + 11 F^V + 6 F^ E 3 + 3 A^EF^ 

-ZA {2)2 B {1) F (0) + 13,4 (2) F^E 2 -6F^ 2 B^ F (0) +3F (i) F (0)2 + AE 2 B {1) F {0) 
_ 7E 2 F (2) A {i) _ 4E 2 F (2) F {o) +8E 2 B (i) A (i) _ 2A {3) E 2 A {1) - A {3) E 2 F {0) + 3 F {2) B {1) F (0)2 

-6 V°> A< J > + 2 EFW A^ 2 + EF™ F<°> 2 - 6 EB™ A< J > 2 + F 4 + 3 F^ B^ F^ A< J > 
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/(24F<°> (aW+fW) (2 A (1 ^ + F^) (3 A^ + F^)) 
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